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Abstract 

A minimal model for mechanical motion of the left ventricle is proposed. The model 
assumes the left ventricle to be a harmonic oscillator with two distinct phases, 
simulating the systolic and diastolic phases, at which both the amplitude and the 
elastic constant of the oscillator are different. Taking into account the pressure 
within the left ventricle, the model shows qualitative agreement with functional 
parameters of the left ventricle. The model allows for a natural explanation of heart 
failure with preserved systolic left ventricular function, also termed diastolic heart 
failure. Specifically, the rise in left ventricular filling pressures following increased 
left- ventricular wall stiffness is attributed to a mechanism aimed at preserving heart 
rate and cardiac output. 
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1 Motivation 



1.1 General 



The left ventricle (LV) is a highly complex mechanical system. Left ventricular 
function is determined by both internal and external factors such as myocar- 
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dial contractility, myocardial stiffness, preload, afterload, rhythm and pace, 
size, shape etc [1]. All these factors combine to regulate left ventricular mo- 
tion, which may be characterized by three basic distinguishing attributes: (1) 
the motion is periodic, including a diastolic expansion phase and a systolic 
contraction phase, (2) the period is divided unequally between the diastolic 
and the systolic phases, the former being generally twice as long as the lat- 
ter [2] and, (3) the forces generated by the left ventricle are in accordance 
with Frank-Starling law of the heart [1], ], which states that the larger the LV 
volume, the greater the energy of its contraction and hence diastolic stretch 
actually increases LV contractility. 

Motivated by the fact that LV dysfunction and heart failure are associated 
with alterations in the mechanical motion of the left ventricle [3,4,5] , sim- 
ulating this motion has become a major challenge. The use of both theoret- 
ical tools (such as finite element geometrical construction[6], fluid-dynamics 
based computation [7,8], Immersed Boundary method [9] etc.) and experimen- 
tal tools (such surgically implanted track markers [10,11], magnetic resonance 
tagging [12,13] etc.) have been utilized for such simulation. As a common 
factor, many of the theories attempt to include an abundance of factors influ- 
encing the left ventricle, to achieve quantitative resemblance between theory 
and experimental data. While some success has been achieved in describing 
normal LV function, the models are far from adequate in describing LV dys- 
function or the effect of cardiac assist devices, and further progress is needed 
before such models may be utilized for clinical applications [14]. 



1.2 Diastolic heart failure 

Heart failure (HF) is a clinical syndrome that can result from structural or 
functional cardiac disorders, impairing the ability of the ventricle to fill with 
or eject blood [15]. The clinical manifestations of HF are dyspnea and fatigue, 
which may limit exercise tolerance, and fluid retention, which may lead to 
pulmonary and peripheral oedema. Diastolic heart failure (DHF) is the clinical 
syndrome of heart failure associated with preserved LV ejection fraction (an 
index of LV systolic function and contractility, defined as the ratio of LV stroke 
volume to end diastolic volume). DHF accounts for approximately 40% of heart 
failure cases [16,17] and carries a significant morbidity, comparable to that of 
systolic heart failure [18,19]. Diastolic abnormalities of the LV may be related 
to increased myocardial stiffness and impaired relaxation [20,21,22], leading 
to an essentially mechanical dysfunction, resulting in the inability of the left 
ventricle to fill with blood at low filling pressures during diastole [23,24]. 

The mechanisms underlying increased myocardial stiffness can be divided into 
factors that are intrinsic or extrinsic to the myocardium [25]. Myocardial fac- 
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tors include cellular factors, such as impaired calcium homeostasis and changes 
in Titin sarcomeric proteins, which act as viscoelastic springs that provide a re- 
coiling force during diastole [26], and extracellular factors, including changes in 
extracellular matrix morphology (i.e. fibrosis) leading to increased myocardial 
stiffness [27]. Relaxation is an energy- dependent process, and abnormalities 
in cellular energy supply and utilization can lead to impaired relaxation [28]. 
Both abnormalities may result in a physiological state in which diastolic LV 
pressures are elevated, leading to elevated left atrial and pulmonary venous 
pressures, exercise intolerance [29,30] and acute pulmonary oedema [31,32]. 

Despite its prevalence and significant morbidity there are currently very few 
models that can simulate DHF as an integral part of the mechanical motion 
of the LV [33]. Such simulation would be beneficial for promoting our under- 
standing of the underlying pathophysiology of DHF and simulating treatment 
modalities, and as such provided much of the motivation for developing a 
minimal model for LV motion. 



1.3 Scope 

In this paper we develop a qualitative minimal approach to describe the motion 
of the LV. Generalizing the simple linear harmonic oscillator, we qualitatively 
mimic the full ventricular cycle, with all three basic distinguishing attributes 
described above. Further more, our model naturally describes the mechanics 
of DHF, and thus provides a baby-step towards a more rigorous model for LV 
function in DHF. 

The rest of this paper is organized as follows. In Sec. 2 we describe and justify 
our model. In Sec. 3 we discuss the relation between our model and DHF, and 
Sec. 4 is devoted to a summary. 



2 The model 

2. 1 Preliminaries 

As a mathematical introduction, in this section we discuss the simple harmonic 
oscillator, describe the general solution and show its relation to the properties 
of the LV. A simple harmonic oscillator (e.g. a mass attached to a spring, an 
elastic band etc.) is described by Hooke's law, that is the larger the displace- 
ment from equilibrium the stronger the restoring force acting on the body. In 
the simplest case, the restoring force is proportional to the displacement, and 
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thus Newton's law is of the form 

r"(t) = -Kr(t) , (1) 

where K is the proportionality (elastic) constant between the force and the 
displacement, its physical meaning is the amount of force required to stretch 
the oscillator by one unit length (per unit mass). 

The solution of Eq. (1) depends on the sign of K. If K is negative, the general 
solution is 

r(t) = r + exp(V^Kt) + r_ exp(- v /Z ^t) , (2) 

where r + and r_ are constants determined by the initial conditions of the 
oscillator. Notice that this solution is "unstable", that is r(t) will increase 
indefinitely with the time t. Thus, this solution is inappropriate for describing 
an oscillating body. 

For a positive K, the solution is periodic, 

r(t) = r cos(v / ^t + 0) , (3) 

where the radial amplitude r and the phase <j) are constants, which are deter- 
mined by the initial conditions of the oscillator. The frequency of oscillations 
is given by = \f~K, and the period (i.e. the duration of a single cycle) is 
t s = 2n/0. We notice that if r(t) stands for the LV radius, than this solution 
already possesses two of the LV characteristics described in the introduction. 
Obviously, the motion is periodic, and the force (given by the right hand side 
of Eq. 1) is in agreement with Starlings law. 

However, the third characteristic of the LV - the unequal division of the cardiac 
cycle - is unfulfilled by the simple harmonic oscillator. Following, we introduce 
a generalization to the simple oscillator, based on a physiologic mechanism, 
which inherently incorporates this feature. 

2.2 Basic model 

In this section we describe our generalization of the simple harmonic oscilla- 
tor. We start by examining the physiologic origin of systolic and diastolic LV 
function. 

The transition from diastole to systole stems from excitation-contraction cou- 
pling. The trigger for systole is the electrical excitation of cardiac muscle cells 
(action potential), which raises the free intra-cellular Ca 2+ concentration, thus 
activating cellular contraction. Contraction leads to elevation of ventricular 
wall tension, elevation of intra-ventricular pressure and eventually ejection of 
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blood from the ventricle, when a positive pressure gradient is achieved. For di- 
astolic relaxation to occur, intra-cellular Ca 2+ concentration must be actively 
reduced, leading to myocardial relaxation [34]. 

In essence, this can be simplified to a model in which systole begins when 
the LV wall is abruptly stiffened (isovolumic contraction) leading to LV pres- 
sure elevation and systolic contraction, after-which the ventricle is suddenly 
relaxed (isovolumic relaxation), followed by a pressure decline and diastolic 
filling. Even a model based on this simplification would require consideration 
of additional factors such as, e.g. environmental effects, blood flow (rheolog- 
ical) effects, spatial form of the LV, energy considerations etc.). This type of 
analysis is probably beyond mathematical rigour. 

We thus suggest a phenomenological model for LV function. We describe the 
LV as a cylindrical elastic membrane. The model thus consists of a harmonic 
oscillator, for which both the elastic constant K (and thus the frequency) and 
the radial amplitude r differ between the systolic and the diastolic phase. 
Since the systolic (diastolic) phase are defined by a decrease (increase) of 
the LV radius r(t), both the elastic constant and the amplitude will depend 
on the sign of the derivative, r'(t) . This will mimic the sudden stiffening 
and relaxation of the cardiac muscle due to changes in intra-cellular Ca 2+ 
concentration. The force induced by the intra- ventricular pressure, F, is given 
(in the cylindrical approximation) by F(t) = 2irr(t)P(t), where P(t) is the 
internal pressure (divided by the LV wall mass and height). Assuming that in 
the systolic and diastolic phases the pressures P s and Pd are approximately 
constant, we obtain the following equation for the radius of the LV, 

Y"(t) + K s (r - r 0)S ) + 2TxrP s = 0, r'(t) < 

(systolic phase) 
r"(t) + K d (r - r 04 ) + 2nrP d = 0, r'(t) > 

(diastolic phase) 

where K d and K s are the effective elastic constants at the diastole and systole 
respectively, and r 0iC z and r , s are the amplitudes at the diastole and systole, 
respectively. Eq. 4 is subjected to the initial conditions 

r(0) = r max ; r'(0) = (5) 

where r max is the maximal radius of the LV, and to the following conditions: 

(1) r(t) is a continuous function. 

(2) The motion is periodic; that is, if T is the period of the function, then 

r(0) = r(T) = r max . 
Let us start by examining each branch of the solution separately, starting with 
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Eq. (4) for the the systolic phase. Demanding that the solution is stable yields 
the condition K s > 2nP s . Denoting J^'p = r s and taking in account the 
initial conditions (Eq. (5)), we have the solution 

r(t) = (r max - r ) cos (^J K s - 2nP s t^ + r . (6) 

The motion is governed by Eq. (6) as long as the first derivative of the solution 
is negative. The derivative changes sign (from negative to positive) at t s = 
VK s -2irP s > for which r (^) = (r max - r ) cosvr + r = 2r - r rnax . Thus, r min = 
2ro — r max (where r m i n is the minimal LV radius) is the initial condition for the 
motion in the diastolic phase, i.e. the solution for the diastolic part of Eq. (4) 
is subjected to the initial conditions, 

r(t s ) = r min ; r'(t s ) = 0. (7) 

Denoting r d = ^^^p d ' ^ ne som tio n fc> r the diastolic stage is given by 

r(t) = (r min - r d ) cos (y/K d - 2nP d (t - t s )^j + r d . (8) 

As in the systolic phase, this branch of the solution is valid when t is in 
the interval t s < t < t s + -j====. Denoting t d = n j\jK d — 2nP d , S = 

\JK S — 2nP s and d = \JK d — 2irP d we obtain the solution for a single cycle, 

r ^ _ ) ^ Tmax _ r «) cos (0st) +r s ; < t < t s ^ 
ij'min f'd) COS (0d(t - t a )) + r d ; t s <t<t s + t d 

From continuity at t — t s and periodicity (i.e. r{t s + t d ) = r(0)) we obtain the 
relation 

Trnax Train / 1 n \ 

r s = r d = . (10) 



We thus obtain a full solution for the motion of the LV, characterized by the 
four parameters t s , t d , r max and r min . The ratio between t d (the duration 
of diastole) and t s (duration of systole) is t d /t s pa 2. The parameters r max 
and r m i n may vary and are determined by both the elastic constants and the 
pressures in the different phases. 

In Fig. 1 we plot the LV volume for several cycles, actual LV data taken 
from a healthy human subject using non-invasive magnetic resonance imaging 
technique, adopted from K. Kindberg [13] (dashed line) and a fit with Eq. (9) 
(solid line). As seen, the model yields fair agreement with the data, and results 
in the ratio between the diastolic and systolic times, t d /t s = 2.152. 
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Fig. 1. Dashed line: LV volume, data taken from a healthy human subject using 
non-invasive magnetic resonance imaging technique, adopted with permission from 
Ref. [13]. Solid line : LV volume, oc irr(t) 2 , as extracted from the model, as a function 
of time (in units of a single cycle period, ibeat)- The fit between the model and data 
yields a ratio between the diastolic and systolic times, td/t s = 2.152. 



3 Diastolic heart failure 



As described in Sec. 1.2, diastolic abnormalities in DHF patients include in- 
creased LV stiffness and elevated LV filling pressures, leading to the symptoms 
of heart failure. In this section we intend to answer two questions: 

(1) Can we simulate the increased LV stiffness within our model? 

(2) Can this model predict and explain the increase of pressure? 

For the first question, the answer seems fairly simple: In the suggested model, 
LV stiffness is represented by the coefficients K s and Thus, in terms of the 
suggested model, the disease can be simulated by increasing by an amount 
AK. 

In order to answer the second question we note that in DHF patients there 
seems to be no significant change in the heart rate, and thus no real change in 
the duration of the diastole (or systole) [32,35]. Since the duration of the phases 
is defined by the frequencies S and 0d, from their definition we have that 0d 
remains constant only if the term Kd — 2nPd remains constant. Thus, if 
is increased by a portion Ai\", then P d must also increase by Hence, this 
model indeed predicts the rise in pressure. In other words, once the diastolic 
stiffness of the LV wall increases, a compensation is achieved by increasing the 
pressure, in order to prevent a significant change in heart rate. 
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4 Summary 



In summary, a naive harmonic model for the mechanical motion of the left 
ventricle was developed. The model consists of a generalized harmonic oscil- 
lator, for which both the radial amplitude and the elastic constant depend on 
the direction of the motion, yielding both the systolic and diastolic phases of 
the cardiac cycle, and allows for a qualitative simulation of left ventricular 
motion. 

Despite its simplicity, the model yields a natural explanation for diastolic 
heart failure, in which the left ventricular wall is stiffened and diastolic filling 
pressures are increased. Our model suggests that the rise in diastolic filling 
pressures is required in order to maintain ventricular volume when the heart 
rate is not significantly changed, as is the clinically observed case [32]. 

One can postulate a correlation between the model and physiological phe- 
nomena as follows: Maintaining cardiac output (CO) is essential for normal 
function. The CO is a product of heart rate and stroke volume. In DHF pa- 
tients, the stroke volume remains unchanged [36,35] and hence slowing heart 
rate should be prevented, or CO will be reduced. However, when our model 
is considered, if LV volumes are to be maintained despite increased stiffening 
(again, as is the clinical finding), LV pressures have to be increased or heart 
rate decreased. It is thus concluded that the increased LV pressures are the re- 
sult of a forced compensation, aimed at maintaining heart rate and CO when 
the LV wall is stiffened. 

Finally, we note that while the model presented here is very much the minimal 
model for LV function, additional phenomenological details, such as ventricular 
pressure- volume loops and atrial parameters may be added in future studies. 
Such additions would render the model analytically unsolvable, yet numerical 
solutions may be possible, and may advance our understanding of different 
physiological phenomena. 
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